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SPECIFICATION 

TITLE OF THE INVENTION 

Channel Estimation method for a Mobile Communication 

5 System 

FIELD OF THE INVENTION 

The present invention particularly relates to an ideal 
channel estimation method for a communication apparatus at a 
10 receiving end, which is used with a wireless communication 
system, 

BACKGROUND OF THE INVENTION 

FIG. 2 illustrates the configurations of wireless 

15 communication apparatuses and wireless communication system to 
which the present invention applies. A wireless communication 
system comprises a wireless communication apparatus at a 
transmitting end 201 and a wireless communication apparatus at 
a receiving end 202 and provides phone calls and data 

20 transmissions via a radio propagation path 203. Although FIG. 
2 and subsequent figures depict only one of two links, such as 
a downlink, the present invention is equally applicable to both 
an uplink and downlink. 
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In the communication apparatus at a transmitting end 201, 
the data to be transmitted is first encoded by a encoder 207 
in a transmitter 205. For wireless communication, 
convolutional coding and turbo coding are frequently used. 
5 Encoded information is entered into a QAM (Quadrature Amplitude 
Modulation) modulator 206 and modulated. A baseband signal 
derived from modulation is converted to a radio- frequency band 
signal by a radio unit 208 and then transmitted. In the wireless 
communication apparatus at a receiving end 202, on the other 

10 hand, a radio signal received from the radio propagation path 
203 via antenna 210 is converted to a baseband signal by a radio 
unit 215. The baseband signal is first demodulated by a QAM 
demodulator 213 in a receiver 211 and then decoded by a decoder 
214. In the QAM demodulator 213, a demodulation process is 

15 performed in accordance with the channel estimation result 
produced by a channel estimator 212. 

As regards encoding for the wireless communication 
apparatus at a transmitting end, it is assumed, as shown in FIG. 
3, that a turbo coding method is employed at a code rate (R) 

20 of 1/3 and at a constraint length (K) of 4 (dotted lines apply 
for trellis termination only). Turbo coding is defined as a 
third-generation mobile communication coding method by 3GPP (3''*' 
Generation Partnership Project) Specification TS 25.212 and one 
of the widely used coding methods. 
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In a QAM modulator of the wireless communication 
apparatus at a transmitting end, long-studied multilevel 
modulation (QAM modulation) is conducted in order to increase 
the transmission efficiency. In the case of 64-QAM, 6 bits can 
5 be transmitted per symbol. As shown in FIG. 4 (401), 6 bits 
(S5, S4, S3, S2, SI, and SO) are grouped into a 3-bit {S2, SI, 
and SO)-dn-phase component (I component) and a 3-bit {S5, S4, 
and S3) quadrature component (Q component). Gray coding is 
conducted so that the difference between adjacent symbols is 

10 1 bit. Transmission takes place with the signal constellation 
shown in FIG. 4 employed* 

In a radio propagation path, the QAM-modulated outgoing 
signal is affected by a combination of amplitude variation and 
phase rotation, which is unique to a radio propagation and called 

15 "fading". To achieve proper demodulation and decoding in the 
wireless communication apparatus at a receiving end, it is 
therefore necessary to estimate such a variation correctly as 
indicated in FIG. 5. It is assumed that the QAM signal whose 
signal constellation is as shown in 501 is transmitted and 

20 subjected to fading having an amplitude variation of G and a 

phase rotation of 9 as shown in 502. The resulting signal 
constellation is as shown in 503. At the receiving end, the 
pilot signal or other signal transmitted for channel estimation 
is used to estimate the amplitude variation G and phase rotation 
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6 in the channel. Further, the complex conjugates of the 
obtained channel estimation result are multiplied together to 

implement -9 rotation and provide phase rotation correction as 
shown in 505. For a QAM or other signal on which information 
5 is superposed in the direction of the amplitude as well, it is 
also necessary to properly estimate the value of amplitude 
variation G for demodulation purposes. 

To achieve channel estimation described above ^ the pilot 
signal for channel estimation is generally prepared separately 

10 from a data signal, which carries information. The pilot signal 
is characterized by the fact that the signal receiving end knows 
what type of signal is transmitted. Therefore, the received 
pilot signal itself indicates the channel situation. If, for 
instance, the pilot signal has a signal constellation of (1,0) 

15 or is expressed as 1 x exp(jO), the received signal is G x exp(j 

8) . It is obvious that the received signal indicates the channel 
condition. When a signal whose quadrature component only is 
subjected to positive/negative inversion is generated from the 

received signal, the resulting signal is G x exp(-j0). It is 
20 therefore easy to obtain a signal that is required for invoking 

(-6) rotation, which is necessary for demodulation. 

There are two pilot signal transmission methods. One 
method provides parallel transmission as shown in FIG. 6 and 
simultaneously transmits a data signal and pilot signal on 



5 



separate conununication channels. The other method, which is 
shown in FIG. 8, uses the same communication channel for a data 
signal and pilot signal while periodically inserting a pilot 
signal. 

5 Channel estimation will now be described in detail while 

explaining how the use of a pilot signal value is timed. 

When a first pilot signal is transmitted in parallel with 
a data signal as shown in FIG. 6, the number of data received 
within the same period of time varies because the data signal 

10 transmission rate is higher than that of the pilot signal. In 
the example shown in FIG. 6, one square corresponds to one signal 
symbol . While one symbol of the pilot signal P3 ( 601 ) is received, 
four symbols (di, dj, dj, andd4) of a data signal (602) are received. 
Therefore, the simultaneously received pilot signal P3 is used 

15 as the channel estimation result for demodulating the four 
symbols (dj, dj, d^, and d4). Thus, the received P3 is used as 
the channel estimation result. 

If the channel is affected by fading, noise, or the like 
at the time of P3 when the channel estimation method shown in 

20 FIG. 6 is used, the communication quality may deteriorate due 
to improper channel estimation. To avoid this problem, a method 
shown in FIG. 7 is used. When this method is used, a plurality 
of pilot signals prevailing before and after an estimation time 
are weighted as shown in 701 and added together (702) for use 
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as the channel estimation result. In this case, however, simple 
weighting coefficients , such as Wl : W2 : W3 : W4 : W5 = 1:2:3:2:1, are 
used so as to merely give greater weighting coefficients to 
signals closer to a data signal (703). 
5 If, on the other hand, pilot signals are periodically 

inserted using the same channel as shown in FIG. 8, the pilot 
signals 801, 803 exist before and after a data signal 802. 
Therefore, the intended purpose can easily be achieved by 
calculating the average value of the pilot signals 801, 802 and 
10 applying the average value, as the channel estimation result, 
to the entire data signal 802 between the pilot signals 801, 
803. 

If the pilot signals are affected by fading, noise, or 
the like during the use of the channel estimation method shown 

15 in FIG. 8, the communication quality may deteriorate due to 
improper channel estimation. To avoid this problem, a method 
shown in FIG. 9 can be used- This method weights a plurality 
of pilot signals prevailing before and after an estimation time 
as shown in 901 and adds them together ( 902 ) for use as the channel 

20 estimation result. In this case, however / simple weighting 
coefficients, such as W1:W2:W3:W4 = 1:2:2:1, are used so as to 
merely give greater weighting coefficients to signals closer 
to a data signal (903). An example of "W1:W2:W3:W4 = 1:2:2:1" 
use is introduced in "Channel Estimation Scheme using the Plural 
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Pilot Blocks for DS-CDMA Mobile Radio" (Ando, et al; Institute 
of Electronics, Information and Communication Engineers Report 
RCS96-72). 

The feature of the channel estimation methods described 
. 5 above is that the channel estimation result derived from the 
information about one or more pilot signals is applied to the 
demodulation of a plurality of data signals (symbols). When 
a conventional channel estimation method is used with a QAM or 
other modulation method by which the amplitude and phase carry 

10 information, highly accurate channel estimation results are 
required for demodulation as well. Therefore, if the same 
channel estimation result is used for a plurality of data signals 
(symbols), the communication quality significantly 
deteriorates due to incorrect channel variation estimation. 

15 Further, even if a conventional method is used to achieve 

channel estimation by weighting pilot signals prevailing at a 
plurality of times, there are no quantitative calculation 
grounds for weighting coefficients. Therefore, if the 
coefficients are improperly set, the communication quality may 

20 seriously deteriorate. 

As the speed of information transmission in a wireless 
communication system increases, an increasing number of systems 
employ the QAM method as a modulation method for the purpose 
of enhancing the frequency use efficiency. Under these 
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circumstances, it is therefore an object of the present 
invention to provide a channel estimation method for estimating 
a channel efficiently and accurately and establishing 
communication at high quality, that is, providing communication 
5 with excellent error rate characteristics, 

SUMMARY OF THE INVENTION 

To achieve the aforementioned object, the present 
invention improves characteristics by calculating the channel 
estimation results variously for all data signal symbols and 
using them for demodulation* To accomplish the object, the 
present invention also pays attention to the frequency 
characteristics of fading in a radio propagation path and 
enhances the accuracy in channel estimation result calculation 
by eliminating the thermal noise that would increase the 
estimation error. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a block diagram illustrating a fourth embodiment 
20 of a channel estimation method according to the present 
invention; 

FIG. 2 illustrates a system configuration according to 
the present invention; 
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FIG. 3 is a block diagram illustrating an encoder for 
a turbo coding method, which is one channel encoding method; 

FIG. 4 illustrates a signal constellation for 64-QAM, 
which is one QAM scheme; 
5 FIG. 5 illustrates how fading and detection changes a 

QAM outgoing signal; 

FIG. 6 illustrates a conventional channel estimation 
method for use in situations where a pilot signal is transmitted 
in parallel with a data signal; 
10 FIG. 7 illustrates a conventional channel estimation 

method that uses pilot signals prevailing at a plurality of times 
in situations where pilot signals are transmitted in parallel 
with data signals; 

FIG. 8 illustrates a conventional channel estimation 
15 method for use in situations where pilot signals and data signals 
are subjected to time-division transmission; 

FIG. 9 illustrates a conventional channel estimation 
method that uses a plurality of pilot signals prevailing over 
a wide range of time in situations where pilot signals and data 
20 signals are subjected to time-division transmission; 

FIG. 10 illustrates a first embodiment of a channel 
estimation method according to the present invention; 

FIG. 11 is a graph showing a frequency spectrum of fading 
during propagation; 
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FIG. 12 illustrates how a low-pass filter is used for 
high-precision fading estimation in a second embodiment of the 
present invention; 

FIG. 13 is a block diagram illustrating a third embodiment 
5 of a channel estimation method according to the present 
invention; 

FIG. 14 illustrates weighting coefficients for use in 
the third and fourth embodiments of the present invention; 

FIG. 15 illustrates the relationship between pilot signal 
10 sampling intervals and channel estimation intervals in the third 
embodiment of the present invention; 

FIG. 16 illustrates the relationship between pilot signal 
sampling intervals and channel estimation intervals in the 
fourth embodiment of the present invention; 
15 FIG. 17 is a graph showing the results of simulations 

that are performed using the fourth embodiment of the present 
invention; and 

FIG. 18 illustrates a fifth embodiment of a channel 
estimation method according to the present invention. 

20 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[First Embodiment] 

To achieve the object described above, the present 
invention uses a pilot signal sampling frequency that is 
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equivalent to the data signal transmission rate (symbol rate), 
and uses a single pilot signal or a plurality of weighted pilot 
signals to achieve channel estimation for each data signal 
symbol. More specifically, the present invention estimates the 
5 data signal's amplitude variation or phase rotation during 
propagation. 

Pilot signals (PO, PI, P2, P3, P4, ••.) 1001 are 
transmitted or sampled at the same transmission rate (symbol 
rate) as for data signals 1005- In FIG. 10, the pilot signals 

10 PO, P4, P8, P12, and P16 (1002), which picked up at 4-symbol 
intervals, are multiplied by weighting coefficients (1003) and 
added together (1004 ) to calculate the channel estimation result 
for a certain symbol dl (1006) within the data signals 1005. 
An alternative is to use only pilot signal P8, which is received 

15 at the same or closest time, for calculating the channel 

estimation result • It is generally believed that a pilot signal 
received at substantially the same time and via the same 
propagation path as for the data signal symbols to be detected 
should be used as the pilot signal for use in channel estimation. 

20 when the pilot signals for a plurality of symbols before and 
after the above-mentioned time are used, high-precision channel 
estimation can be achieved. While the timing relationship 
between one data symbol and the pilot signal used for channel 
estimation of the data symbol is fixed, the channel estimation 
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result for data symbol d2 (1010), which is next to data symbol 
dl (1006), is calculated by similarly multiplying, weighting 
coefficients (1008) using pilot signals PI, P5, P9, P13, and 
P17 (1007) and adding them up (1009). 
5 When the above channel estimation method is used, the 

channel estimation result can be calculated variously for all 
data signal symbols* As a result, the communication quality 
characteristics (error rate characteristics) can be improved. 

10 [Second Embodiment] 

In a second embodiment, a received pilot signal is used 
as an input and allowed to pass through a low-pass filter (LPF) 
for the purpose of enhancing the channel estimation accuracy. 
The second embodiment will now be described in detail. 

15 In a radio propagation path, fading occurs in which the 

amplitude and phase drastically change (a sort of noise) . When 
this phenomenon is expressed in terms of a frequency region, 
no frequency components exceed the maximum Doppler frequency 
(= fd). The maximum Doppler frequency is proportional to the 

20 carrier wave frequency and mobile station's moving speed. If, 
for instance, the carrier wave frequency is 2 GHz and the mobile 
station's movement speed is 120 km/h, the maximum Doppler 
frequency (f^) is substantially equal to 240 Hz (equivalent to 



13 

240 bps) so that the substantial spectral density is as indicated 
by 1101 in FIG. 11. 

Meanwhile^ the data transmission rate of, for instance, 
the third-generation cellular phone is 384 kbps. The 
5 transmission rate used in an actual radio propagation path is 
one order of magnitude higher. In an example shown in FIG. 11, 
the data transmission rate is set at 289.5 kHz for simulation 
purposes (1103). In this case, however, a rate of 289.5 kbps 
is more than 1200 times higher than an f^ value of 240 Hz. 

10 During propagation, however, thermal noise (white noise) , 

which is proportional to the bandwidth of a radio propagation 
path, is applied in addition to the above-mentioned fading as 
indicated by 1201 in FIG. 12. Conventional channel estimation 
is conducted while a broadband noise component 1201 is added 

15 to fading-induced amplitude variation and phase rotation ( 1202 ) . 
However, the fading does not have a component having a frequency 
higher than the maximum Doppler frequency as described earlier. 
Therefore, when a pilot signal is allowed to pass a low-pass 
filter having a passband for passing components having a 

20 frequency lower than the maximum Doppler frequency as indicated 
by 1104, the fading signal passes through as is, but the majority 
of thermal noise within the remaining band is eliminated. It 
is needless to say that the thermal noise existing within the 
fading frequency region remains without being eliminated ( 1203 ) . 
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However, when compared to the entire thermal noise initially 
applied during propagation, the amount of the remaining thermal 
noise can be considered extremely small because the data 
transmission rate is high and the fading is one order of magnitude 
5 later (the above filter can reduce the amount of thermal noise 
to about 1/600 in the above example in which the value fd is 
240 Hz and the transmission rate is 289.5 kbps ) . When a low-pass 
filter appropriate for the maximum Doppler frequency of a fading 
invoked during propagation is applied to a pilot signal used 
10 for channel estimation as described above, it is possible to 
eliminate unnecessary noise and achieve channel estimation with 
high accuracy. 

[Third Embodiment] 

15 The description of a third embodiment relates to the 

configuration of a low-pass filter that is used in the second 
embodiment. For calculating the channel estimation result, 
which varies from one data signal symbol to another as mentioned 
in the description of the first embodiment, the pilot signal 

20 transmission or sampling intervals are the same as those for 
the data signal. In the first embodiment, pilot signals at a 
plurality of times are weighted and added together. However, 
the simplest method for achieving this object is the use of an 
FIR filter 1301 shown in FIG. 13. The FIR filter comprises a 
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plurality of delay devices 1303, multipliers 1304 for 
multiplying weighting coefficients for respective data, and an 
adder 1305 for adding up the obtained multiplication results. 
As shown in FIG. 14, various FIR filter coefficients are set 
5 as weighting coefficients in accordance with an impulse response 
1402, which is determined by a Fourier transform of desired 
low-pass characteristics 1104 (1404). It can be said that the 
FIR 1301* is used to form a low-pass filter having frequency 
characteristics 1104* One feature of the present invention is 

10 that the weighting coefficients are set so as to provide low-pass 
filter characteristics appropriate for channel variations. 
When attention is given to the weighting coefficient values, 
another feature of the present invention is that negative 
weighting coefficient values may be used as indicated by 1404 

15 in FIG. 14. 

[Fourth Embodiment] 

A fourth embodiment uses the Nyquist Theorem to perform 
channel estimation result calculations with high accuracy. In 
20 the fourth embodiment, the channel estimation accuracy attained 
when the pilot signal sampling rate is lower than the data 
transmission rate (data signal symbol rate) is equal to the 
accuracy attained when the pilot signal sampling rate is the 
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same as the data transmission rate. The fourth embodiment will 
now be described in detail. 

According to the Nyquist Theorem, the original signal 
can be reproduced when the sampling frequency is more than two 
5 times the maximum frequency of the original signal. When 
sampling is conducted at 480 Hz (equivalent to 4 80 bps ) or higher 
in the case of a fading in which the value f^ is 240 Hz 1102 
as shown in FIG. 11, the fading can be theoretically reproduced. 
The case indicated in FIG. 11 satisfies the Nyquist Theorem's 

10 conditions because the data transmission rate is 289.5 kbps 
(1103) . The third embodiment described earlier also satisfies 
the Nyquist Theorem when the pilot signal sampling intervals 
are rendered equal to the data transmission rate, as indicated 
in FIG. 15. The pilot signal 1501 is sampled at a sampling 

15 frequency equivalent to the data transmission rate 1502, and 
an FIR filter 1301 is operated to pass a frequency region 
appropriate for fading in order to calculate a channel 
estimation result 1503 at intervals equivalent to the data 
transmission rate. 

20 When the maximum Doppler frequency f^ for fading is 240 

Hz, the operation of the third embodiment permits signal 
reproduction when sampling is conducted at a rate of 480 Hz or 
higher. However, the data transmission rate is as high as 289.5 
kbps. In this case, the pilot signal sampling rate can be 
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reduced to about two times the maximum Doppler frequency without 
impairing the communication quality. 

As shown in FIG. 16, therefore, the sampling rate is 
reduced, for instance, to 1/4 the data transmission rate, that 
5 is, sampling is conducted every four data symbols, and the value 
"0" is inserted at the remaining three sampling times. The 
resulting signal 1601 is subjected to interpolation and then 
entered into an FIR filter that is used in the third embodiment. 
The data transmission rate is adequately high relative to the 

10 maximum Doppler frequency of fading. Therefore, even when the 
pilot signal sampling rate is reduced to 1/4 the data 
transmission rate, the Nyquist Theorem is satisfied. Further, 
since the employed FIR filter is the same as used in the third 
embodiment, all the frequencies contained in the fading pass 

15 through. Consequently, the channel estimation result 1602, 
which is output as a result of the above operation, is 
substantially equal to the channel estimation result 1503 shown 
in FIG. 15. The employed channel estimation method is shown 
in FIG. 1. The pilot signal 1001 is sampled (1002) at a certain 

20 sampling frequency (which is, however, at least two times the 
maximum Doppler frequency of fading) and used as an input signal 
for channel estimation. Next, a switch 102 is operated so as 
to generate a signal 1601, into which the value "0" is inserted 
in accordance with the channel estimation result output 
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intervals {= data transmission rate), and then an FIR filter 
provided with weighting coefficients having low-pass . 
characteristics indicated by 1404 is operated. The FIR filter 
comprises a plurality of delay devices 103, multipliers 104 for 
5 performing multiplications with respective weighting 

coefficients, and an adder 105 for adding up the obtained 
multiplication results. The output 106 of the FIR filter 
operating at a data symbol rate becomes a channel estimation 
result for a data signal symbol at a specific time. Subsequently, 

10 a detector (107) uses the channel estimation result to achieve 
data signal detection (phase rotation correction, etc.), and 
a demodulator 108 demodulates the data signal that is modulated. 

The error rate characteristics prevailing when the fourth 
embodiment of a channel estimation method were evaluated by 

15 means of computer simulation. For simulation purposes, turbo 
encoding was conducted at a code rate (R) of 1/3 and at a 
constraint length (K) of 4 in FIG. 3, and 16-QAM multilevel 
modulation was conducted to perform a transmission at 20 kbps 
(the prevailing QAM symbol transmission rate was 16 K symbols /s) . 

20 Further, the pilot signal was transmitted at 1600 Hz (1.6 kbps) . 
It was assumed that the pilot signal power was 10 times higher 
than the data signal power. For channel estimation purposes, 
the pilot signal was sampled every ten data signal symbols, and 
nine Os were inserted for agreement with the data transmission 
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rate. A 100-tap FIR filter was used as a low-pass filter having 
a passband of frequencies up to 240 Hz. FIG. 17 shows the error 
rate characteristics prevailing when f^ was 240 Hz. 

The solid line represents an error rate characteristic 
5 boundary value that can be attained when ideal channel 

estimation is accomplished without error. If there is residual 
noise 1203^ which cannot be eliminated by the FIR filter, the 
error rate characteristics deteriorate due to noise-induced 
error in channel estimation. The resulting deteriorated 

10 characteristics, which are indicated by the dotted line, may 
be interpreted as a limit that prevails when the low-pass 
characteristics appropriate for fading characteristics are 
implemented by the FIR filter. The line with triangular marks 
indicates the characteristics that prevail during the use of 

15 a conventional channel estimation method disclosed by NTT Docomo 
in 1996. FIG. 17 reveals that the use of the present embodiment 
of an FIR filter having low^pass characteristics appropriate 
for fading characteristics introduces a characteristics 
improvement of about 0,2 dB over the conventional method. 

20 

[Fifth Embodiment] 

When a multilevel modulation is conducted in a first 
embodiment, the employed channel estimation result varies from 
one data signal symbol to another in order to ensure high 



20 



efficiency and accuracy in channel estimation. In the fifth 
embodiment, however, the channel estimation result is updated 
every several symbols so that the channel estimation method used 
in the first to fourth embodiments applies at intervals of 
5 several symbols. The channel estimation result intrinsically 
varies from one data symbol to another. Therefore, the 
characteristics deteriorate to a certain extent when the same 
channel estimation result is applied to several symbols. 
However, the frequency of channel estimation result calculation 

10 can be reduced. Further, when an FIR filter is used, the number 
of taps can be reduced from the viewpoint of hardware 
configuration (equivalent to a decrease in the number of 
multipliers). As a result, increased simplicity is provided 
as an advantage. One example is shown in FIG. 18. The example 

15 presented in FIG. 18 assumes that six taps are provided for an 
FIR filter 1801 for use in updating the channel estimation result 
for each symbol, and that the distribution of employed weighting 
coefficients is as indicated by 1802. Further, when it is 
assumed that the same channel estimation result may be used for 

20 two adjacent data symbols, the use of three multipliers suffices 
because the two adjacent symbols use the same weighting 
coefficient for multiplication. The resulting distribution is 
as indicated by 1804 (the dotted line in 1804 represents 
weighting coefficients that are equivalent to the 



21 ' 

characteristics prevailing when the number of taps is 3 and 
implemented when the number of taps is 6 as is the case with 
1801). Although a certain degree of characteristics 
deterioration results from the use of different weighting 
5 coefficients, the number of multipliers in this example can also 
be reduced to half. 

The pilot signal for use with the present invention may 
be either a common pilot signal for transmission from a 
transmitting communication apparatus to a plurality of 

10 receiving communication apparatuses or a specific pilot signal 
for transmission from a transmitting communication apparatus 
to a specific receiving communication apparatus * If pilot signal 
sampling is conducted at a rate lower than the data signal symbol 
rate as is the case with the fourth embodiment, the employed 

15 pilot signal may also be a pilot signal that is on the same channel 
as for data signals and periodically inserted between data 
signals. 

Even when multilevel-modulated data is transmitted at 
a high speed, the present invention makes it possible to achieve 
20 channel estimation with high efficiency and high accuracy, 
thereby improving the communication quality characteristics. 



